Two-speed phase dynamics in Si(lll) (7x7)-(lxl) phase transition 
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We propose a natural two-speed model for the phase dynamics of Si(lll) 7x7 phase transition 
to high temperature unreconstructed phase. We formulate the phase dynamics by using phase-field 
method and adaptive mesh refinement. Our simulated results show that a 7x7 island decays with 
its shape kept unchanged, and its area decay rate is shown to be a constant increasing with its 
initial area. LEEM experiments concerned are explained, which confirms that the dimer chains and 
corner holes are broken first in the transition, and then the stacking fault is remedied slowly. This 
phase-field method is a reliable approach to phase dynamics of surface phase transitions. 

PACS numbers: 68.35.-p, 05.10.-a, 68.37.-d, and 05.70.-a 



Introduction It is well known that bulk terminated 
lxl semiconductor surfaces usually are unstable against 
reconstruction at low temperatures and some of the 
reconstructed surface phases transit to the lxl ones 
at high temperatures Si surfaces are well stud- 
ied because silicon plays the centering role in modern 
computer industry 0, 0, H[ and potential silicon-based 
spintronicsQ. The most important Si surface phase 
is the Si (111) 7x7 reconstructed surface [||. It is the 
ground state structure of Si(lll) surface and transits 
to an unreconstructed lxl structure above T c =1125 
K0, 1 i 1 El E i EI EE 0- In situ 
scanning-tunneling-microscopy (STM) observations 
supported by Monte Carlo simulations [l3l 17 1, indicated 
that the stacking faulted triangle unit is the basic build- 
ing block in forming the 7x7 surface phase. The phase 
transition is believed to be of first order 11, IE 13. 14. 
15, 12, IS ! IE], although earlier data implied a continuous 
phase transition [3]. It is of much interest to clarify the 
phase dynamics during the phase transition. Low-energy 
electron microscopy (LEEM) imaging is a powerful ap- 
proach to experimentally measure time-dependent sur- 
face phases [3, HE Hl|. Recent LEEM results showed 
that big 7x7 islands always decay faster than small 
ones, with area decay rates increasing with their initial 
areas It seams like that the area decays have some- 
thing like momentum. This phenomenon is very intrigu- 
ing. It means that the phase transition dynamics and the 
microscopic processes meanwhile are very much complex, 
what we have known about them may be just a tip of the 
iceberg, and therefore much more investigation is in need. 

Here we investigate the time-dependent phase dynam- 
ics during the phase transition of Si(lll) 7x7 islands to 
the unreconstructed high temperature phase. Through 
analyzing experimental results and the atomic structures 
of both surface phases, we propose a two-speed model 
for the phase dynamics. We formulate the phase dynam- 
ics using phase field method, famous for various growth 
and solidification issuesfH, 0, HJ EE HE HE HE], and 
adaptive mesh refinement technique [291]. Our simulated 



results indicate that the dimer chains and corner holes 
are broken first in the phase transition, and then follows 
the slower remedying of the stacking fault. Our simu- 
lated images and linear area decays of 7 x 7 islands are in 
agreement with the experimental results, which supports 
the two-speed model. This is a reliable approach to un- 
derstand the phase dynamics of the Si(lll) 7x7 phase 
transition and other surface phase transitions. 

Experimental clues The bulk terminated Si(lll) lxl 
surface and the dimers-adatoms-stacking-fault (DAS) 
model of Si (111) 7x7 reconstructed surface phase are 
shown in Fig. 1 [H, @] . This lxl surface is not experi- 
mentally realized. The real experimental phase at high- 
temperature is a "1 x 1" phase that is formed by cover- 
ing the bulk terminated lxl surface with 0.25 monolayer 
(ML) of fast-moving adatoms [IE Hit The 7x7^"lxl" 
phase transition needs extra 0.17 ML of adatoms because 
the former has 0.08 (4/49) ML more adatoms than the 
bulk terminated lxl surface. On the experimental side, 
the appearing of the 7x7 LEED pattern, or the seventh- 
order spots, means that the "lxl" -^7x7 transition hap- 
pens, and on the other hand the absence of it means 
that the 7x7 structure is destroyed [32j|. The adatoms 
stay at T4 sites above the second- layer atoms in the 7x7 
surface, and some of them moves to H3 sites above the 
forth-layer atoms when the surface transits to the "lxl" 
phase [33J. The movement of adatoms from T4 sites to 
H3 sites or vice versa is easy because it can be realized 
without breaking any bond. The formation of the dimer 
chains, corner holes, and stacking fault is considered es- 
sential to the 7x7 reconstruction [34| . The presence of 
the LEEM image means that the key 7x7 factors, the 
dimer chains and corner holes, still exist. Without the 
key factors, the 7x7 phase is destroyed. The resultant 
region without dimer chains and corner holes, although 
still having stacking- fault, cannot be distinguished from 
the "lxl" phase by means of LEEM experiment, being 
'ghost' region to LEEM imaging |35j|. 

Model and phase-field realization Therefore, we believe 
that the atoms reorganization during the 7x7^ "lxl" 
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FIG. 1: (color online). The atomic structures of bulk termi- 
nated Si(lll) lxl surface (left) and the DAS model of Si(lll) 
7x7 reconstructed surface (right). The top four layers are 
shown for both structures, and in addition, the 12 adatoms 
are included for the latter. The 7x7 unit is divided into the 
faulted (F, upper) half and the unfaulted (U, lower) half by 
the dash line. 



transition can be described by two basic processes: (a) 
the destroying of the dimer chains and corner holes, and 
(b) the remedying of the stacking- fault. The latter pro- 
cess should take more time because it requires collective 
movement of many atoms concerned. We adopt phase 



field method £2, 23 



24 l25l |26| ] to simulate the phase 
dynamics during the 7x7^ "lxl" transition. We use 
phase-field variables cf> and £ together to describe the 
7x7 phase with respect to the "lxl" phase. </> is used 
for the aspect of the dimers and corner holes according 
to process (a) and £ for the aspect of the stacking fault 
according to process (b). Both <j> and £ are functions of 
the time t and two-dimensional space coordinates (x,y), 
and are made to have two stable values, 1 and -L as is 
done in usual phase-field simulations 22, 2^, 24, 25|, l26j |. 
The adatoms are described with another variable u to be 
defined in the following. Therefore, the complete 7x7 re- 
gion, with not only the dimer chains and corner holes but 
also the stacking- fault, is described by 0=1 and £=1; and 
the complete "lxl" region by <j>=-l and £=-1. <f*=l, re- 
gardless of £, reflects the presence of the dimer chains and 
corner holes, the key 7x7 character detected in LEEM 
imaging experiment. £=1 and <f>=-\ together describe the 
temporary 'ghost' regions during the 7x7 -"lxl" phase 
transition. The governing equations can be expressed as 



T0 ~dt = w ^ 2 ^ ~ ~ x 3<t>{4>)u - X' ' g<t>(4>)g(Q 
Te§ = ^ 2 V 2 £-/ e (£)-A' 5 (^(£) 



(1) 



where V 2 is defined as + Jp-, and the functions f v (v) 
and g-qiv) {il=4>, £) are the derivatives of functions f(rj) 
and g{rj), respectively. We define f(rj) as exp[(r/ 2 — l) 2 ], 
rather than usual {rf — l) 2 , because the A' terms are 
fifth-order polynomials of <fi and £. This choice not only 
keeps the two stable values, +1 and -1, and the parabolic 
shape in their neighborhoods without adding extra ex- 
treme values, but also guarantees the stability of the sys- 
tem against the A' terms. The and rj describe the 
characteristic evolution times of the phase fields. The 
V 2 term and the time derivative one together determines 
the evolution rate of the phase-field, with D v =w 2 /t v act- 
ing as controlling parameter. In addition, w 2 determines 
the transition zone of the phase-fields in the asymptotic 
regime. Here we have < because process (b) is 
slower than process (a). The A term describes the in- 
teraction between u and <j>, reflecting the growth or the 
shrinkage of LEEM-detectable islands with the help of 
adatoms. The A' term describes enhanced evolution of 
<j>=l (£=1) islands in the presence of £=1 (0=1). We use 
,g (77) =(2 + 377 — ?y 3 )/3 because the growth and shrinkage of 
the 0=1 region takes place only by the u exchange along 
its boundary and is enhanced only with the existence of 
£=1 regions. 

Computational methods and parameters The variable 
u describes the local density difference of adatoms be- 
tween the "lxl" structure and the 7x7 one. Because 
the adatoms move very quickly, we suppose that u does 
not directly depend on space positions, but is a func- 
tion of the phase field 0, u=uq[1 — tanh(0)/ tanh(l)J/2. 
This means that u is zero deeply in the 0=1 islands [l9(, 
equal to Mo deeply in the <f>=-l regions, and in between in 
transition zones where takes values between 1 and -1. 
We take uq— 1.33nm -2 in terms of experimental adatom 
measurements P, |3l|. Other parameters are taken in 
terms of basic time interval dt=0.08s and basic grid 
length dir=1.77 nm. The latter is determined in terms 
of the area of the minimal 7x7 equilateral triangle. Be- 
cause the half 7x7 unit cell is an equilateral triangle, we 
conserve the threefold anisotropy by using w,j,=wih(9$), 
W^=W2h(6^), T^—r\h 2 (Q^), and r^—Tih 2 (9^). The func- 
tion h(9 v ) (ry=0, £) is defined as h(9 v )=l + ecos(39 v ), 
where e parameters the anisotropy, and the directional 
function 9, q is defined as arcsin(— r/y/ \Jvx + Vy) ~ tt/2, 
where the r] x and rj y denote the derivatives of rj with re- 
spect to x and y. We take e=0.3, w\=W2=dx 1 T\=dt, and 
T2=3dt in our simulated results. The ratio D^/ D^=Q.33 
is the relative rate of process (b) with respect to process 
(a). The interaction constants A and A' are 2.34g?x 2 and 
2.225, respectively. All our simulated results are robust 
enough against changes of the parameters. 

We use an adaptive mesh refinement method [2^] to 
perform effectively our phase-field simulations in two di- 
mensional space. The adaptive mesh refinement enables 
us to simulate a large spatial scale (~ /jm) within an ac- 
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ceptable computational time, and to get enough details 
about the regions where the phase fields change drasti- 
cally without adding too much computational time. 
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FIG. 3: (color online). Simulated decay behaviors of five </>=l 
islands with initial areas 50=0.1210, 0.0764, 0.0431, 0.0211, 
and 0.0093 /im 2 (from top to bottom) within the single phase- 
field model (left) and the two-phase-field model (right). In the 
left panel, the circle indicates the simulated island areas and 
the curve is fitted with the function So exp(— at— (3t 3 ' 5 ). In the 
right panel, the filled (hollow) circles indicate the simulated 
areas of 0=1 (£=1) islands and the lines are linearly fitted 
0=1 results. The inset shows the f=l areas over the whole 
simulation time. 




FIG. 2: (color online). Simulated images of the and £ fields 
within the two-phase-field model. White triangles denote the 
0=1 islands and black background the <f>=-l region in (a), (b), 
(e), and (f), and red (grey) triangles denote the £=1 islands 
and blue (black) background the £=-1 region in (c), (d), (g), 
and (h). The time interval is 200, 400, or 300 time steps 
between (a) and (b), (b) and (e), or (e) and (f), respectively, 
where 100 time steps are equal to 8 seconds [ly]. 

Main findings through simulation The simulated im- 
ages are presented in Fig. [2l The whole time sequence of 
the 7x7-^ "lxl" phase transition is shown in the eight 
panels (a-h), among which (a), (b), (e), and (f) for 
and (c), (d), (g), and (h) for £. At the beginning, the 
£=1 region is always enclosed in a larger 0=1 region, as 
shown in the four panels (a)-(d). The £=1 region must 
be right in size with respect to the 0=1 region in order 
to achieve the linear area decay behavior of the 7x7 re- 
gions. Too small £=1 region is not enough to drive the 
decay from the exponential behavior, and on the other 
hand too large £=1 region overdrives the area decay from 
the linear behavior. The right area difference satisfies the 
condition that S® — S9 is a linear function of 5°. This 
can be understood by considering the fact that when the 
sample is quenched, the slow faulted stacking lags the 
formation of the dimcrs and corner holes. 

The area decays for five islands are presented in the 
right panel of Fig. [3l The filled circles are the areas 
of 0=1 islands and the hollow ones those of £=1 re- 
gions. The area decays are linear with time and the de- 
cay rates are 0.0060, 0.0046, 0.0031, 0.0021 and 0.0014 



(/im 2 per 100 time steps) for the initial areas: 0.1210, 
0.0764, 0.0431, 0.0211, and 0.0093 ^m 2 , respectively. It 
can be concluded that the area decay rate of a 7x7 is- 
land increases with its initial area, being independent of 
its current area. An initially larger island always decays 
faster than an initially smaller one even when their sizes 
become equivalent to each other. The area decay rate of a 
7x7 island remains the same until the island disappears. 

If using only one phase-field variable, we get area decay 
results shown in the left panel of Fig. [3] The filled circles 
show the areas of five islands with the same five initial 
areas as those of the two-phase-field model. They can 
be fitted with the function So exp(— at — /3t 3 ' 5 ), with the 
parameters a and [3 decreasing with So- Area decay rates 
of different 7x7 islands increase with their current areas, 
not directly depending on their initial areas. The rate 
is 0.0047 [ive? per 100 time steps at S=0.04/zm 2 (upper 
dash line), and reduces to 0.0024 /Ltm 2 per 100 time steps 
at 5=0.01/im 2 (lower dash line). This behavior is not 
compatible with the experimental LEEM results [ijj], and 
therefore the two phase-field variables are necessary to 
explaining the experiment. 

Comparison with experiment In Fig. 2]we compare our 
simulated area decay results (0=1) in terms of the two- 
phase-field model with experimental data concerned. In 
the inset we present the relation of the area decay rate 
(R) with the initial area (So). We obtain a linear func- 
tion, R = 0.00015 + 0.00515 . With this function, the 
area decay rates are obtained by comparing the simulated 
island areas with experimental ones, or vice versa. The 
experimental island area S can be described with a linear 
function of the time t, S = So — Rt. Because initial island 
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FIG. 4: The two-phase-field simulated area decay (lines) of 
five 7x7 islands compared with experimental LEEM data 
(dots)pjj. Inset: the simulated decay rate (R in units of 
10~ 4 fxm 2 1 s) vs the initial area So (fJ.m 2 ), where the squares, 
derived from Fig. 3, are fitted with a linear function. 

area So can be determined more accurately in LEEM 
experiment, we obtain So by fitting experimental area 
data. The fitted S are 0.0976, 0.0737, 0.0475, 0.0216, 
and 0.0090/im 2 , and the corresponding area decay rates 
are 6.4, 5.2, 3.9, 2.6, and 2. Ox lO^V" 2 /*- The theo- 
retical linear behaviors are in good agreement with the 
experimental data[19(. Therefore, our theoretical model 
is very good for understanding the experimental results. 

Conclusion We have proposed a natural two-speed 
model for the phase dynamics of the phase transition 
of Si(lll) 7x7 islands to unreconstructed high tempera- 
ture "lxl" phase. We formulate the phase dynamics by 
using phase-field method and adaptive mesh refinement 
technique. Our simulated results show that a 7x7 island 
decays through step flow and the triangular shape is kept 
all the time. The decay rate of the island area is shown 
to be approximately a constant increasing with the initial 
area only. The LEEM experiments are explained quan- 
titatively. This in return supports our conclusion: the 
corner holes and dimer chains are broken first and the 
remedying of the stacking- fault is slower and takes longer 
time. Therefore, the phase dynamics is elucidated. This 
phase-field theory is a reliable approach to studying the 
phase dynamics of general surface phase transitions. 



This work is supported by Nature Science Founda- 
tion of China (Grant Nos. 10774180, 90406010, and 
60621091), by Chinese Department of Science and Tech- 
nology (Grant No. 2005CB623602), and by the Chinese 
Academy of Sciences (Grant No. KJCX2.YW.W09-5). 



* Electronic address: bgliu@mail.iphy.ac.cn 
[1] W. Monch, Semiconductor Surfaces and Interfaces, 

Springer, Berlin 2001. 
[2] E. Kasper and D. J. Paul, Silicon quantum integrated 

circuits, Springer, Berlin 2005. 
[3] J. B. Hannon et al, Nature 440, 69 (2006). 
[4] Y. Sugimoto et al, Nature 446, 64 (2007). 
[5] I. Appelbaum et al, Nature 447, 295 (2007). 
[6] K. Takayanagi et al, Surf. Sci. 164, 367 (1985). 
[7] T. Hoshino et al, Phys. Rev. B 51, 14594 (1995). 
[8] T. Hoshino et al, Phys. Rev. Lett. 75, 2372 (1995). 
[9] W. Shimada and H. Tochihara, Surf. Sci. 526, 219 (2003). 
[10] F. J. Giessibl et al, Science 289, 422 (2000). 
[11] J. Kanamori and Y. Sakamoto, Surf. Sci. 242, 102 and 
119 (1991). 

[12] S. Hasegawa et al, Phys. Rev. B 47, 9903 (1993). 
[13] T. Kato et al, Surf. Sci. 416, 112 (1998). 
[14] C.-W. Hu et al, Surf. Sci. 487, 191 (2001). 
[15] H. Hibino et al, Phys. Rev. B 72, 245424 (2005). 
[16] P. A. Bennett and M. W. Webb, Surf. Sci. 104, 74 (1981). 
[17] M. Itoh, Phys. Rev. B 54, 5873 (1996); 56, 3583 (1997). 
[18] W. Telieps and E. Bauer, Surf. Sci. 162, 163 (1985). 
[19] J. B. Hannon et al, Nature 405, 552 (2000). 
[20] J. B. Hannon et al, Phys. Rev. Lett. 86, 4871 (2001). 
[21] J. B. Hannon, J. Tersoff, and R. M. Tromp, Science 295, 
299 (2002). 

[22] A. Karma and W. J. Rappel, Phys. Rev. E 53, R3017 

(1996); Phys. Rev. Lett. 77, 4050 (1996); Phys. Rev. E 

57, 4323 (1998). 
[23] A. Karma, Phys. Rev. Lett. 87, 115701 (2001). 
[24] F. Liu and H. Metiu, Phys. Rev. E. 49, 2601 (1994). 
[25] O. Pierre-Louis, Phys. Rev. E. 68, 021604 (2003). 
[26] Y.-M. Yu and B.-G. Liu, Phys. Rev. E 69, 021601 (2004); 

Phys. Rev. B 70, 205414 (2004); Phys. Rev. B 73, 035416 

(2006). 

[27] D. D. Vvedensky, J. Phys. CM 16, R1537 (2004). 

[28] J. W. Evans, P. A. Thiel, and M. C. Bartelt, Surf. Sci. 

Rep. 61, 1 (2006). 
[29] N. Provatas, N. Goldenfeld, and J. Dantzig, Phys. Rev. 

Lett. 80, 3308 (1998); J. Comp. Phys. 148, 265 (1999). 
[30] Y.-N. Yang and E. D. Williams, Phys. Rev. Lett. 72, 1862 

(1994). 

[31] Y. Fukaya and Y. Shigeta, Phys. Rev. Lett. 85, 5150 
(2000). 

[32] W. Telieps and E. Bauer, Ultramicroscopy 17, 57 (1985). 
[33] S. Kohmoto and A. Ichimiya, Surf. Sci. 223, 400 (1989); 

T. Suzuki et al, Phys. Rev. B 59, 12305 (1999). 
[34] T. Hoshino et al, Surf. Sci. 394, 119 (1997). 
[35] E. Bauer, Rep. Prog. Phys. 57, 895 (1994). 



